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We present a method for scheduling observations in small field-of-view tran- 
sient targeted surveys. The method is based in maximizing the probability 
of detection of transient events of a given type and age since occurrence; it 
requires knowledge of the time since the last observation for every observed 
field, the expected light curve of the event and the expected rate of events in 
the fields where the search is performed. 

In order to test this scheduling strategy we use a modified version of the 
genetic scheduler developed for the telescope control system RTS2. In par- 
ticular, we present example schedules designed for a future 50 cm telescope 
that will expand the capabilities of the CHASE survey, which aims to detect 
young supernova events in nearby galaxies. We also include a brief description 
of the telescope and the status of the project, which is expected to enter a 
commissioning phase in 2010. 



1 Introduction 

With a new generation of observatories dedicated to studying the time do- 
main in astronomy [1, 2, 3, 4 , our understanding of astrophysical transient 
phenomena will be significantly improved. The diversity of known families of 
transient events will be better understood thanks to improved sample sizes 
and better data, and new types of transient events will be likely discovered. 

These observatories will include large field-of-view, large aperture tele- 
scopes, which will scan the sky in a relatively orderly fashion, but also net- 
works of small field-of-view, small aperture robotic telescopes that will scan 
smaller areas of the sky in a less predictable way. 

The smaller robotic telescopes are ideal for studying very short-lived tran- 
sients, e.g. gamma ray bursts (GRBs), but also to do detailed follow up stud- 
ies of longer lived galactic (e.g. cataclysmic variables, planetary systems) and 
extragalactic (e.g. supernovae) transient events. Moreover, they constitute a 
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relatively inexpensive tool to obtain reduced cadences, of the order of days, 
in relatively small areas of the sky which are of special interest, e.g. nearby 
galaxies. 

Here, we present a scheduling strategy that maximizes the probability 
of finding specific types of transient phenomena, or the expected number of 
events, at different times since occurrence. In Section [2] we derive the proba- 
bility of finding one or more of these events, as well as the expected number 
of events. In Section |3] we show the results obtained with this method and 
discuss its implications. Finally, in Section|3]we give an overview of the future 
50 cm telescope that will expand the capabilities of the CHASE survey and 
which will use the scheduling method presented in this work. 

2 Detection probabilities of transient events 

This discussion will be limited to well known types of events in targets with 
known distances. We assume that the light curves of every transient event is 
composed of a monotonically increasing early component, followed by a mono- 
tonically decreasing late component. We will show how to compute detection 
probabilities for individual difference observations, as well as for sequences 
of observations to pre-defined targets. With this information, we will discuss 
how to build observational plans that maximize the detection of events with 
certain characteristics. 

Expected numbers vs probabilities 

The probability of having exactly k occurrences of an event, in a time interval 
where A occurrences are expected, is: 



which means that the probability of zero occurrences of the event is: 



which we would like to minimize. Note that for small values of A, the prob- 
ability of detecting at least one event should be a better indicator of a good 
schedule than the total expected number, but since I — e~ A m 1 — (1 — A) = A, 
in practice this can be ignored. For big values of A the total expected number 
of events should be most of the time a better indicator of a good schedule 
than the probability of finding at least one event. For the purpose of this 
discussion we will use probabilities, but it is easy to change the formulation 
to the expected number of events, as we will show later. 



P(k,X) 




(1) 



P(0,A) = e-\ 



(2) 



Scheduling in targeted transient surveys and a new telescope for CHASE. 



3 



Detection probabilities of individual events 

Let us assume that the events remain detectable for a time r and that their 
rate of occurrence is R. Consider also the case when we look at a target twice 
to generate a difference image, with a time interval or cadence, At. 

If each event remains visible for r years, we would like to know what is 
the interval where an event which was not seen in the first observation could 
occur and be detected in the second observation. 

Let us also assume that the event was not seen in the first observation, 
performed at time t\ , and that we make a second observation with a cadence 
At, i.e. at time ti = t\ + At. Defining At! as the minimum between At and 
t, then the time interval where new transients can occur and be detectable in 
the second observation will span from £2 — At' and ti . This is because short- 
lived transients only have a time r to remain visible, which could be smaller 
than At. Hence, the expected number of new events that can be detected will 
be the rate of occurrence times the former time interval. Using equation ((2|), 
the probability of no events occurring in this interval and no detections being 
made, Pnd, will be: 

Pnd = cxp{-i? mm(At, r)}. (3) 

Let us now assume that the event can only be detected a time to after its 
occurrence, that it remains visible for a time r and that we are only interested 
in events younger than r age (see Figure [T]). 



T age 








T ' \ 




Time 


T - t 



''age l D 



Fig. 1. The different variables used in this calculation, to is defined as the time when 
the luminosity is above a certain threshold value which makes the object detectable, 
r is defined as the time period where the transient can be detected and r ago is an 
arbitrary age since explosion, with r age — to being the time period where events 
younger than r agc can be detected. 

The event will be detectable younger than r ago only if r agc > to. If this is 
the case, the time period where events not seen in the first observation could 
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occur and be detected in the second observation will now be the minimum 
between At, r ago — to and r (see Figured]). 

Thus, the probability of no events occurring in this time interval and no 
detections being made, P^d > wm ^ e: 



nd — ex P 



-R mm{At, max(r agc — to, 0), t} 



(4) 



With this information, the probability of detecting one or more events in 

,age 
ND • 



the second observation will be simply 1 — Pj asf 
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Fig. 2. Example of how the detection probabilities are calculated. At is the cadence, 
t is the time period when the transient remains detectable and r agc — to is the time 
period when the transient remains detectable while it is younger than T age . In this 
case, the cadence is shorter than r, which means that some transient events of the 
type of interest could have already been seen in the first observation, but r agc — to is 
shorter than the cadence, which means that all events younger than r ago could not 
have been seen in the first observation. The minimum between these three quantities 
should be multiplied by the rate of events in order to compute the expected number 
of events and compute the detection probabilities. 



Cadence choice 

Using the formula above, we could try maximizing the probability of detection. 
For a fixed target, this can only be done decreasing the cadence, At, as long 
as the number of targets that are observed in the sample is not compromised 
significantly. 

It is easy to see that, if r ago > to, increasing At from zero to larger values 
will increase the probability of detection only while At < min(r agc — to,r). 
For bigger cadences the probability will stay constant. Thus, a natural choice 
for the cadence would be At = min(r agc —to,r). 



Scheduling in targeted transient surveys and a new telescope for CHASE. 



5 



If larger cadences were chosen, the probability of detecting events younger 
than T ag0 would not change, but not all detected events would be guaranteed 
to be younger than r ag0 , which could be a problem if a fast age estimation is 
required. On the other hand, larger cadences would increase the probability of 
detecting older events, up to when At > r, when the probability of detecting 
an event of any age would remain constant too. 

Choosing At = min(T ag0 — to, t) has the added benefit that if T age is smaller 
than the rise time and the object type was known, the absolute magnitude 
of the event could be used as an age estimator. This is because the event 
would be guaranteed to be rising at detection time and the magnitude-age 
relation would be single-valued. In reality, there could be events of different 
types simultaneously occurring which would make the age determination only 
useful in a statistical sense. 

In general, At will be a function of the distance and light curve of the 
variable object to be detected, but also of the desired age of detection, T ago . 
Hence, for young and bright objects with extended light curves, the cadence 
should be set to at least r ag0 — to if we want to increase the cadence while 
maximizing the detection of events with a given age. 

However, it is not always easy to repeat the observations with a fixed 
cadence. Bad weather, the change of position of the targets throughout the 
year or the appearance of other objects of interest, among many reasons, may 
cause the cadence between observations to vary. 

An alternative strategy is to let the cadence adapt individually in a se- 
quence of observations in order to maximize the detection probabilities. 

Detection probabilities for a sequence of observations 

Now, we compute the probability of not detecting any new events in a sequence 
of observations, P^d Total- 

We note that for no events to be detected, each individual observation 
must result in negative detections, i.e. we have: 

P nd 6 , Total = n i exp[-ff min{At\ max(r agc - t , 0), t 1 } 

= exp -Si FC mm{At\ max(r agc - t , 0), t 1 } , (5) 

where the indices indicate different targets. Thus, the probability of detecting 
one or more new events will be: 

P D S Totai = 1 - cxp[-i7 t FT mm{At\ max(r agc - t , 0), t 4 }] . (6) 

With this formula, we recover the expected number of events in the en- 
tire observational sequence, which is the term inside the exponential, i.e. in 

^D, Total — 1 - ex P( _ ^D S Total)' 

*D%tai - R l ram{At\ max(r agc - ij , 0), t 1 } (7) 



6 Francisco Forster, Nicolas Lopez, Jose Maza, Petr Kubanek, G. Pignata 

is the expected number of new detected events with the required age. Thus, 
we can use either equation ([5]) or ([7]) to determine the fitness of individual 
schedules, but we recommend using equation (J7J- 

Limited number of targets 

It is possible that the number of targets available for detecting new events 
with a given age is too small, i.e. assuming a fixed exposure time and cadence 
for all observations, that the number of visible targets where r age > t is 
smaller than the length of the night divided by the exposure time. 

For very short-lived transient surveys this is not a problem, since even 
with relatively small cadences At > r, and the probability of detection in an 
individual observation would be p w R r, i.e. it would be almost independent 
of the cadence or how many times we observe a target per night. 

In relatively long-lived transient surveys, i.e. time-scales of days or longer, 
we would not want to repeat targets in a given night. This is because when 
At < £ agc — to the probability of detection in an individual observation is 
P w R At. Thus, many observations to a given target in a given night would 
be almost equivalent to observing the target once per day or once every few 
days in terms of probabilities, but with a significantly higher cost on the 
resources and preventing the telescope from observing other targets. 

In general, the number of targets for a given cadence should be of the order 
of the fraction of time that we want to spend in that sample per night, /, times 
the total number of observations per night, N cxp , times the cadence, At. The 
detection rate would be approximately the multiplication of this number with 
the typical rate of occurrence, R. 

Hence, a possible strategy would be to order targets by the time that it 
takes for the events of interest to be detectable, to, depending on their distance 
and extinction, select a detection age according to scientific criteria, and then 
group the targets according to the resulting cadence and sample sizes. This is 
summarized in the following table: 

Table 1. For a desired detection age (T age ) we show a possible choice of cadence 
in days (At), the maximum sample size consistent with this cadence and the ap- 
proximate detection rate if the suggested cadence and maximum sample sizes were 
used, to corresponds to the time for the transient events to become visible since 
its occurrence, A^xp is the number of observations per night and R is the rate of 
events expected in each field. We assume that in all cases the time that a transient 
remains visible, r, is bigger than the cadence chosen. For the opposite case, At can 
be replaced by r in the last two columns below: 



Age 


Reference Cadence 


Sample size 


Approx. detection rate 


T~agc 


At = (r ag c - t ) 


/ 7V cxp At 


/ Nexp At (R) 



Scheduling in targeted transient surveys and a new telescope for CHASE. 7 
2.1 Genetic algorithm 

As discussed above, one can let the cadence vary from observation to observa- 
tion and from object to object. For an ideal schedule, we would like to select 
the optimal combination of cadences that can adapt to unexpected changes 
of the observational plan. For this, we use the probability of detection, or the 
expected number of events, of a sequence of observations as the fitness indi- 
cator and we use a genetic algorithm to find the best available observational 
plan for the following night. This can reflect unexpected changes to the ob- 
servational plan in a daily basis, and can be extended to fractions of a night 
optimizations if necessary. 

We have used the genetic algorithm implemented in RTS2 5 , taking into 
account the cadence to each target (At) and the distance, event rate, height 
above the horizon and sky brightness, all of these reflected in the quantities 
to and r, to build the observational plan. 

The distance between targets is also taken into account indirectly. If it is 
too big, the number of visited targets per night, or the number of terms in 
Equation ([7]), will be reduced and the probability of detection will decrease 
accordingly. Similarly, when the targets are too distant, or too close to the 
horizon, or the sky too bright, to will increase and t will decrease, decreasing 
the detection probabilities too. The bigger the event rate in every target, 
the bigger the detection probability, which will favour those targets with the 
biggest intrinsic rates. Finally, the time since last observation will determine 
the cadence, changing the detection probabilities as well. 

In these calculations, the time between targets is computed using the max- 
imum between the slew time and the readout time, which effectively defines 
a disk around each target where the time penalty is constant. Reaching the 
outer circumference of this disk would take exactly the readout time assum- 
ing that the CCD can read out electrons while simultaneously slewing in the 
most efficient trajectory. This is regularly accomplished by RTS2, since it op- 
timizes observations by reading out electrons and moving to a new position 
simultaneously. 

For instance, a readout speed of about 2 sec and a slew speed of 5 deg 
sec -1 define a disk around the previous target of about 10 deg in the sky where 
the time penalty for new targets is the same. In most telescopes the slewing 
movement is accomplished with two independent motors, which makes the 
size and shape of this disk really depend on the initial configuration of the 
telescope before slewing, and whether an equatorial or altazimuthal mount is 
used. 

The details of the genetic optimizer, based on the NSGAII algorithm [TJ, 
are described in detail in [8] . It is worth mentioning that the genetic algorithm 
can handle multiple objectives, which can be used to find the Pareto front of 
optimal values instead of a single solution, e.g. look for multiple detection 
ages, which we have also implemented (see Figure E}. 
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The Pareto front is the locus of solutions in a multi-objective optimization 
problem where one objective cannot be improved without compromising the 
other objective functions. For example, in an optimization problem with two 
objective functions, for every value of one of the two objective functions there 
is an optimal value for the remaining objective function, i.e. the Pareto front 
can be composed by infinite solutions. 



2.2 Calculation of t and r 



In the previous sections we did not include the calculation of the time for an 
event to become detectable, to, and the time that an event remains detectable, 
r. These terms can be computed from empirical light curves of the particular 
event to be detected, and can be stored as functions of the critical luminosity 
above which the object can be detected. 

Thus, the problem is reduced to computing the flux above which the object 
can be detected. To do this, we solve the signal to noise equation for an 
arbitrary value above which we define an object to be detected, e.g. S/N = 5. 
This equation is: 

S/N(t) = 7te ( ) m , (8) 

[7^ D (t)T + 7s ^D Tnpix + 7 2 N npix] V 2 

where S/N{t) is the signal to noise ratio as a function of time, 7te D M are 
the photons per unit time detected by the CCD from the transient event as a 
function of time, T is the exposure time, 7^y D are the photons per unit time 
coming from the sky and detected in one pixel of the CCD, n p ; x is the number 
of pixels used to do photometry and 7rn is the readout noise per pixel of the 
CCD. In general, n p i x is a function of the seeing at the zenith and the angle 
from the zenith. 

Solving the previous quadratic equation for 7te° (*) with a given value of 
S/N and choosing the positive root gives the following result: 



v ,ccd (4.\ _ S/N' 



7te (*) 



2 T 



4 Tipix r 2 CCD T"ri i 

-gjj^YlBN +7sky ±\) 



1/2' 



(9) 



Now, we can include the effect of distance, collecting area, spectral shape and 
CCD characteristics in the following equation: 

7?e D W - ^2 / 7teW V, dv (10) 

where A is the collecting area of the telescope, D is the distance to the object, 
7te(*) i s the number of photons per unit time per unit solid angle per unit 
frequency of the transient event as a function of time since occurrence, r/i, is 
the efficiency with which the photons are captured as a function of frequency, 
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which depends on the reflecting surfaces, intervening lenses, CCD quantum 
efficiency and filters. 

We can write a similar equation for the photons coming from the sky in 
every pixel of the CCD: 

where Ail is the solid angle of one pixel of the CCD, z is the angle from the 
zenith and 7g KY is now the number of photons coming from the sky per unit 
time per unit area per unit solid angle per unit frequency. 

Thus, if we compute 7te D (*)> assuming the object is at a distance D from 
the observer, 7tez)Wj an d 7sky assuming the object is at a given angle 
z from the zenith, 7sky ( z )> f° r a gi ven sky brightness and for a particular 
telescope configuration, we can simply scale the results as follows: 

ZV 2 



=7$£ W [jf) ( 12 ) 

Thus, we can now compute the times when the object becomes detectable and 
when it is no longer detectable, to and to + t: 



to, t + ■ 



7tE,D 




4 77 ■ 7 C t CD (0) \ 1/2 
"17rn H — — 1 i 



S/N 2 



(14) 



where ^7te£) J ^ s ^ ne Averse of the function computed in equation (|12p. 

which should have two solutions for a transient which is composed by an early 
monotonically increasing component followed by a monotonically decreasing 
late component. Importantly, the inversion of 7xez)o mu st be performed only 
once, and can be stored numerically in a table, e.g. in logarithmic intervals of 
photons per unit time. 

Thus, for a given signal to noise ratio (S/N), which we arbitrarily define 
as the value that gives a detection, a given distance from the source (D), 
exposure time (T), sky brightness (c.f. 7sky ) : seeing (c.f. n np i x ), readout 
noise per pixel (trn) and angle from the zenith (z), we can compute to and 
r, which are necessary for the calculation of the expected number of events 
and the probabilities of detection in an individual target and a sequence of 
observations. 

It is important to note that the detection of objects is sometimes performed 
using individual pixels, in which case we can set n p j x to one, and multiply the 
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term 7te W m equation © by the fraction of photons that fall in the central 
pixel in the position of the object, depending on the seeing conditions, which 
would result in the following modified equation: 



to, to 



„, CCD 
7te,_d 




S/N 2 



{;: 



RN 



7^(0) 
cosz 



T} 



1/2 



(15) 



where / is the fraction of the light from a point source that would fall in one 
pixel in the position of the object, generally a function of the seeing at the 
vertical, the angle from the zenith and the frequency of the photons to be 
detected. 



3 Results and discussion 

In Figures [3] and [4] we present example implementations of the scheduling 
strategy presented in this work with the genetic algorithm used in RTS2. 

Figure [3] shows the probability of finding supernova with a reference 50 cm 
telescope in an observational plan composed of 60 sec individual exposures, 
with simulated cadences and supernova rates in each field. We can see the 
probability increasing with each generation of observational plans and then 
staying constant. Each generation is formed by a population of 1,000 different 
observational plans and the initial iteration consisted of a series of randomly 
generated targets for each observational plan of the population, which were 
crossed and mutated to obtain the best observational plans. 

Figure @] shows the space of optimal solutions when two objective functions 
are used. This is, the Pareto front of non-dominated solutions, or the space 
of solutions where one variable is at its optimal value without compromising 
the other variables. In this simulation we use the objective functions: (1) 
probability of finding supernovae before maximum and (2) the probability 
of finding supernova younger than three days from explosion, using similar 
parameters to those used in the simulation shown in Figure [3] 

Interestingly, we have used the already implemented genetic scheduler from 
RTS2 to find the schedule that maximizes the average height above the horizon 
for our list of targets, or that minimizes the typical distance between targets. 
For both cases, we have found that the probability of detection of the resulting 
schedule is smaller by more than a factor of two with respect to our method, 
which suggests that our strategy is significantly better for finding transient 
objects. 

Thus, the implemented scheduling strategy based on maximizing the prob- 
ability of finding new transient events is able to obtain significantly higher 
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detection probabilities than alternative methods. We were able to build ob- 
servational plans for every night to maximize the probability of detecting par- 
ticular events, or similarly, the expected number of detections. These plans 
were based on pre-defined samples of targets that have characteristic cadence 
and exposure times, and that can easily adapt to unforeseen changes in the 
scheduled observations. 

In order to compute the observational plans with the highest detection 
probabilities, we used the genetic algorithm implemented in the telescope 
control system RTS2, where a multi-objective algorithm selects the optimal 
sequence of observations for our purposes. 

We expect to be able to extend this work to scheduling of coordinated 
networks of robotic telescopes looking for specific types of transient events, or 
looking for many different phenomena if multi-objective optimization is used. 
We also expect to release the implementation in a future version of RTS2 
(http://rts2.org). 

An important question is whether this method is able to recompute the 
optimal observational plan when unexpected changes in the sequence of ob- 
servations occur. In a single computer, with the current implementation of the 
code we cannot think of simple ways of achieving this, since it normally takes 
many hours to find the optimal observational plan or set of Pareto-optimal 
plans in a single PC. However, with faster computers, pre-calculating detec- 
tion probabilities for every target at every time in the night, and given that 
genetic algorithms can be relatively easily parallelised, we expect this to be 
feasible in the near-future. 

Alternatively, one could switch from using optimized observational plans to 
computing the detection probabilities for every available target and choose the 
one with the highest detection probability every time the telescope has finished 
integrating, taking into account the slew and readout time by subtracting the 
expected cost of slewing in term of detection probabilities per unit time for 
the corresponding slewing times. 

Finally, it should be noted that this method is not exclusive for supernova 
transients, but to any transient with well characterized light-curves and with 
well understood target fields. 

4 Application to the new 50 cm robotic telescope for 
CHASE 

The CHASE survey [6] is the most prolific nearby supernova search in the 
southern hemisphere. It finds more than 70% of the nearby (z < 0.3) su- 
pernova in the southern hemisphere, with discovery ages much younger than 
competing surveys (see Figure [5]). CHASE uses a fraction of the time available 
in four of the six PROMPT telescopes pQ located in CTIO. 

In order to expand the capabilities of CHASE and to have a better control 
over the scheduling of the observations, we are in the process of purchasing 
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Fig. 3. Evolution of the probability of detecting new supernova with the number 
of generations. In this example, every generation consists of a population of 1,000 
observational plans, where each observational plan contains hundreds of 60 sec ex- 
posures to different targets. For the simulation, we have used the gold sample of 
galaxies of the CHASE survey [6], where the distances were computed from reces- 
sion velocities. Event rates were created randomly at the beginning of the simulation 
to reproduce the typical rates expected for supernova explosions. The cadences were 
randomly generated to reproduce the typical values expected in the CHASE survey. 

and installing a 50 cm robotic telescope that will join the other PROMPT 
telescopes for the SN survey and follow up. 

The telescope will be a 50 cm automated telescope: composed of an optical 
tube, a CCD camera with a set of filters, a mount, a meteorological station, a 
dome and computers for controlling and analyzing the data. It will be located 
in CTIO and remotely controlled from Cerro Calan (Santiago, Chile). It will 
observe hundreds of targets every night with the aim of doubling the observing 
capabilities of the CHASE survey and to try new observing strategies with 
new associated scientific goals. 

The optical tube of the telescope will be a 50 cm aperture Ritchey- 
Chretien design, with a focal ratio of 12, in an open-truss carbon fiber tube 
purchased from the Italian company Astrotech. The camera will be a 2kx2k 
pixels Finger Lakes Proline camera, with a back illuminated, UV enhanced, 
95% peak quantum efficiency Fairchild 3041 CCD. The pixel size will be 0.52" 
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Fig. 4. Locus of Pareto-optimal solutions, or Pareto front, using two objective func- 
tions: the probability of finding SNe before maximum (abscissa) and the probability 
of finding SNe no later than three days after the explosion (ordinate). In this ex- 
ample we have evolved 10,000 generations of 1,000 observational plans each, similar 
to the simulation shown in Figure [3l but with two objectives instead of only one. 
Absolute numbers should not be trusted, since the simulated cadences are too big 
for the values expected in a survey looking for SNe as young as three days after 
explosion. Individual points represent individual observational plans, and the set 
of points are only an approximation to the Pareto front. Once the Pareto front is 
computed, one observational plan from the set of solutions can be chosen according 
to arbitrarily defined criteria. 



and the field-of-view will be 17.6' in side. The camera will be equipped with 
a 12-slot filter wheel with the filters u'g'r'i', Johnson B and V and WFCAM 
Z, purchased from Asahi-Spectra (see transmission curves in Figure [6]). 

The camera was chosen to avoid the potential presence of residual images in 
the imaging of targets, which currently dominate our SN candidate lists with 
the PROMPT telescopes, to obtain a relatively big field-of-view, which would 
allow us to image enough reference stars to do an accurate image alignment 
and subtraction, but also to obtain the best available quantum efficiency, 
which is a cost-effective way of collecting more photons per target. 

The mount will be the Astro-Physics 3600GTO "El Capitan" model, which 
is a German equatorial mount with sub-arcmin pointing errors, and a slew 



14 Francisco Forster, Nicolas Lopez, Jose Maza, Petr Kubanek, G. Pignata 




-15 -10 -5 5 10 15 5 10 15 20 25 30 



Fig. 5. CHASE survey results. Top: the distribution of supernova with declination 
in bins of equal solid angle. Red bins correspond to all nearby SNe found in 2008, 
and green bins correspond to SNe discovered by CHASE in 2008. Bottom-left: the 
distribution of discovery ages of Type la SNe in CHASE. About half of the SNe were 
discovered before maximum. Bottom-right: the distribution of discovery ages of Type 
II SNe in CHASE. The median discovery age was about 5 days after explosion. (G. 
Pignata, private communication). 

speed of about 5 deg sec -1 . The dome of the telescope will be built in Chile 
and is currently in the design phase. 

The scheduling of the observations will be done with the strategy presented 
in this work, and we expect to start collaborations with other groups using 
this scheduler in an integrated fashion. For more information please contact 
the authors. 
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Fig. 6. Approximate transmission curves for the filters implemented in our future 
50 cm telescope. "Open" corresponds to no filters, and includes the effects of atmo- 
spheric extinction, mirror reflectivity and CCD quantum efficiency. The supernova 
search will be likely performed in "Open" mode or with a clear filter, but for the 
follow-up program we will use the filters shown here. It is possible that in the future 
we will include additional filters. 
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